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We analyze baryon number, strangeness and electric charge fluctuations as well as their correla- 
tions in QCD at high temperature. We present results obtained from lattice calculations performed 
with an improved staggered fermion action (p4-action) at two values of the lattice cut-off with almost 
physical up and down quark masses and a physical value for the strange quark mass. We compare 
these results, with an ideal quark gas at high temperature and a hadron resonance gas model at low 
temperature. We find that fluctuations and correlations are well described by the former already 
for temperatures about 1.5 times the transition temperature. At low temperature qualitative fea- 
tures of the lattice results are quite well described by a hadron resonance gas model. Higher order 
cumulants, which become increasingly sensitive to the light pions, however show deviations from a 
resonance gas in the vicinity of the transition temperature. 

PACS numbers: 11.15.Ha, ll.10.Wx, 12.38Gc, 12.38.Mh 



> 

O 
O 



oo 
O 

> 

X 



I. INTRODUCTION 



Fluctuations of conserved charges, like baryon number, 
electric charge and strangeness, are generally considered 
to be sensitive indicators for the structure of (subsets of) 
a thermal medium produced in heavy ion collisions [l[ . In 
fact, if at non- vanishing baryon number a critical point 
exists in the QCD phase diagram, this will be signaled by 
divergent fluctuations of e.g. the baryon number density 
0- 

Under conditions met in current experiments at RHIC 
as well as in the upcoming heavy ion experiments at the 
LHC the net baryon number is small and QCD at van- 
ishing chemical potential provides a good approximation. 
In this region the transition from the low temperature 
hadronic to the high temperature plasma regime is con- 
tinuous and fluctuations are not expected to lead to any 
singular behavior. Nonetheless, they provide direct in- 
sight into the structure of the thermal medium, the rele- 
vant degrees of freedom and their correlations. Further- 
more, enhanced fluctuations provide hints for nearby sin- 
gularities in the QCD phase diagram related to the chiral 
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limit at vanishing net baryon number as well as for a pos- 
sible critical point at physical values of the quark masses 
at non- vanishing net baryon number density Q. 

Away from criticality, i.e. under conditions met at 
RHIC and LHC, indications for the existence of such crit- 
ical points can only show up in higher order derivatives of 
the QCD partition function with respect to temperature 
or chemical potentials. Through the analysis of fluctua- 
tions of conserved charges as well as their higher moments 
and correlations we thus gain insight into the relevant de- 
grees of freedom of the system under consideration and 
at the same time gather information on possible nearby 
singularities in the QCD phase diagram. 

From lattice calculations at vanishing chemical poten- 
tial it is well known that baryon number and strangeness 
susceptibilities are sensitive indicators for the transi- 
tion from the low temperature hadronic regime to the 
high temperature quark gluon plasma. Also in the case 
of a continuous crossover transition, as is the case for 
QCD with physical quark masses, the susceptibilities rise 
rapidly in the transition region 0, H, S 0, Hi- This 
shows that on the scale of the temperature the (quasi)- 
particles carrying the quantum numbers under consider- 
ation (B, S) are heavy at low and light at high temper- 
ature. 

In calculations with two light, dynamical quark de- 
grees of freedom (2-flavor QCD) it also has been shown 
that baryon number and electric charge fluctuations in- 
crease and that their fourth moments start to show pro- 
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nounced peaks in the transition region from low to high 
temperature In fact, higher order cumulants of e.g. 

baryon number fluctuations, become increasingly sensi- 
tive to the singular behavior in the vicinity of the chiral 
phase transition at zero mass and vanishing chemical po- 
tential. Starting from the 6th order cumulant baryon 
number fluctuations will diverge in the chiral limit and 
are expected to reflect critical behavior in accordance 
with the 0(4) symmetric universality of the chiral tran- 
sition. 

It also has been noted that the analysis of correlations 
between different quantum numbers or flavor channels 
will provide insight into the quasi-particle structure and 
the relevant degrees of freedom in QCD at high temper- 
ature [l(J. In lattice calculations this has been analyzed 
in quenched [TO] and 2-fiavor QCD 0,10. Recent 
lattice calculations for the QCD equation of state per- 
formed with dynamical light and strange quark degrees of 
freedom [3, EH now also allow to perform a systematic 
analysis of these correlations among different quantum 
number channels including effects arising from the treat- 
ment of strange quarks as dynamical degrees of freedom 

dm. 

We present here results from lattice calculations of 
baryon number, electric charge and strangeness fluctu- 
ations in QCD with dynamical light and strange quark 
degrees of freedom 1 . The results are based on calcula- 
tions with an improved staggered fermion action (p4- 
action) that strongly reduces lattice cut-off effects in bulk 
thermodynamics at high temperature. The values of the 
quark masses used in this calculation are almost physical; 
the strange quark mass, m s , is fixed to its physical value 
while the light up and down quark masses are taken to 
be degenerate and equal to m s /10. This is about twice 
as large as the average up and down quark masses re- 
alized in nature. We obtained results from calculations 
performed with two different values of the lattice cut-off, 
corresponding to lattices with temporal extent N T = 4 
and 6. This allows us to judge the magnitude of sys- 
tematic effects arising from discretization errors in our 
improved action calculations. The spatial volume has 
been chosen to be V 1 I 3 T = 4, which insures that finite 
volume effects are small. 



tials can be expressed in terms of chemical potentials 
for baryon number (hb), strangeness (us) an d electric 
charge (jiq), 
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Moments of charge fluctuations, 5Nx = Nx — (Nx), 
with X = B, Q or S and their correlations are then 
obtained from derivatives of the logarithm of the QCD 
partition function, i.e. the pressure, 



P _ 1 



T 4 VT 3 
evaluated at ^b.q.s = 



In Z(V,T, /j, s , n Q ) , 
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with fix = fix/T. While the first derivatives, i.e. baryon 
number, electric charge and strangeness densities, vanish 
for fiB,Q,s = 0, their moments and correlation functions 
with i + j + k even are non-zero. The basic quantities 
we will analyze here are the quadratic and quartic charge 
fluctuations 2 , 
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and the correlations among two conserved charges, 



All 



f 



VT 3 



(NxN Y ) 



(4) 



(5) 



These quantities have been evaluated in the temper- 
ature interval 0.8<T/T C <2.5 on lattices of size 16 3 x 4 
and 24 3 x 6, respectively. On the 16 3 x 4 lattice we also 
calculated 6th order expansion coefficients, 



II. FLUCTUATIONS AND CORRELATIONS; 
COMPUTATIONAL SETUP 



At vanishing baryon number (B), electric charge (Q) 
and strangeness (S) fluctuations of these quantities can 
be obtained by starting from the QCD partition func- 
tion with non-zero light and strange quark chemical po- 
tentials, pL u ,d,s = fJ-u,d,s/T. The quark chemical poten- 
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3Q(N X ) 3 ) 



(6) 



The gauge field configurations, that have been used to 
evaluate the above observables, had been generated pre- 
viously in calculations of the QCD equation of state [lij 
and the transition temperature [1 61 ] . In these calculations 
the strange quark mass has been tuned close to its physi- 
cal value and the light quark masses have been chosen to 
be one tenth of the strange quark mass. This corresponds 



1 Preliminary results of this calculation had been presented at Lat- 
tice 2007 and 2008 [3. 



2 As all expectation values have been evaluated at vanishing chem- 
ical potential, we have SNx = Nx- 
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302 
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363 


350 
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345 


10 


50 



TABLE I: The data samples analyzed on lattices of temporal 
extent N T = 4 and 6, respectively. The columns give from left 
to right the temperature values, the number of configurations 
used in this calculation, the number of trajectories by which 
these configurations are separated and the number of random 
vectors used for the evaluation of traces. 



to a line of constant physics on which the kaon mass is 
close to its physical value and the lightest pseudo-scalar 
mass 3 is about 220 MeV, i.e. the light quark masses used 
in these calculations are about a factor 2 larger than their 
physical values. Further details on the improved gauge 
and staggered fermion actions used in these calculations 
are given in [hH . [l6| . The number of gauge field config- 
urations analyzed varies from about 300 at high to 1500 
at low temperatures. Subsequent configurations are sep- 
arated by 10 to 60 trajectories. The various operators 
contributing to xfjk nave been calculated using unbi- 
ased random estimators 18]. While at high temperature 
already 50 random sources per configuration were suffi- 
cient to get reliable estimates for these observables, we 
used up to 480 random sources below the transition tem- 
perature. Autocorrelations of the operators contributing 
to the quadratic and quartic fluctuations turned out to 
be about 100 trajectories close to T c and to drop quickly 
away from it. Errors on these fluctuations have been de- 
termined by means of the jack-knife procedure. Some 



In calculations with staggered fermions flavor symmetry is bro- 
ken at non-vanishing values of the lattice spacing a. As a con- 
sequence only one of the pseudo-scalar mesons has a light mass 
that is proportional to •J'mi and vanishes in the chiral limit at 
fixed a > 0. Full chiral symmetry with the correct Goldstonc 
pion multiplet is recovered only for a — -» 0. For a study of the 
remaining flavor symmetry violations with the p4-action in the 
quenched case see fl7ll . For dynamical calculations discussed 
here we have calculated one of the non-Goldstone pion masses 
at cut-off values corresponding to the transition region of our 
jV T = 4 and 6 calculations. This gives 700 MeV and 550 MeV, 
respectively. 



details on our data sample are given in Table Q] 

In the following we will present our results using a 
physical temperature scale in MeV. As explained in detail 
in Ref. [l4[ , this scale has been obtained through detailed 
studies of the zero temperature static quark potential 
along the line of constant physics used also for the finite 
temperature calculations. This yields unambiguously the 
temperature in units of the lattice scale tq that charac- 
terizes the shape of the potential at distance Tq, i.e. r$ is 
defined through the relation (r 2 dVqq(r)/dr) r _ ro = 1.65. 
Although ro is not directly accessible in experiments, its 
value is quite well known through comparison with deter- 
minations of, e.g. the level splitting in the bottomonium 
system as well as calculations of light meson decay con- 
stants [HI, [2(| • For the conversion to physical units we 
have used the value r = 0.469 fm (20j . 

In calculations with the improved staggered action the 
transition temperature has been determined on lattices 
with temporal extent N T = 4 and 6 for the values of quark 
masses used here [H, This gave T c = 204(2) MeV 
for N T = 4 and 196(3) MeV for N T = 6, respectively. 
The temperature values given in Tables I-III have been 
obtained from a fit to ro/a. The fitting function is given 
in Eq. 22 of Ref. [14] . We therefore do not quote errors on 
individual temperature values in the table. As discussed 
in Ref. [l4j, these errors are generally below 1%. 



III. FLUCTUATIONS OF LIGHT AND 
STRANGE QUARK NUMBERS 

Before entering a discussion of fluctuations of B, Q and 
S it is instructive to look into fluctuations of the partonic 
degrees of freedom, the light and strange quarks. It has 
been noticed before that close to the transition temper- 
ature fluctuations of the heavier strange quarks are sup- 
pressed relative to those of the light u or d quarks [HJ ■ At 
higher temperatures, however, they approach each other. 
In these earlier calculations, however, the strange quark 
sector has not been incorporated as a dynamical degree 
of freedom in the numerical calculations. Nonetheless, 
the basic observation also holds true in QCD calculations 
with dynamical strange quarks. 

The difference in the behavior of light and strange 
quark fluctuations is apparent from Fig.[TJ where we show 
the quadratic fluctuations of light and strange quark 
numbers (upper figure) and the ratio of these fluctuations 
(lower figure). The quadratic fluctuations rapidly ap- 
proach the Stefan-Boltzmann, ideal gas value, X2 = 1j 
and decrease strongly in the transition region. The cal- 
culations on lattices with temporal extent N T = 4 and 6 
show some cut-off dependence, which is larger for light 
quarks than for strange quarks. Nonetheless, it is appar- 
ent that the approach to the high temperature limit is 
slower for the heavier strange quarks than for the light 
up or down quarks. This is highlighted in the lower part 
of Fig. [1] where we show the ratio of strange to light 
quark fluctuations. This ratio shows less cut-off depen- 
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T[MeV] 
176.0 
186.4 
190.8 
195.4 
202.4 
204.8 
209.6 
218.2 
253.7 
305.2 
444.4 



X2 

0.1656(28) 
0.3530(41) 
0.4298(50) 
0.5144(74) 
0.7095(58) 
0.7689(66) 
0.9452(52) 
1.0040(27) 
1.0556(28) 
1.0543(21) 
1.0158( 5) 



T[MeV] 



173.8 
179.6 
185.6 
195.0 
201.3 
204.6 
211.1 
224.3 
237.7 
278.4 
362.9 
415.8 



X:l 



0.2077(77) 
0.2769(85) 
0.3950(81) 
0.6011(116) 
0.7556(80) 
0.8093(86) 
0.8714(65) 
0.9313(61) 
0.9695(33) 
0.9975(48) 
1.0181(35) 
1.0163(18) 



xj_ 

0.0922( 9 
0.2186(20 
0.2660(30 
0.3190(46 
0.4464(38 
0.4918(49 
0.6577(54 
0.7371(25 
0.9078(25 
0.9851(18 
0.9913( 5 



Xi 



0.1167(35 
0.1643(32 
0.2341(45 
0.3743(93 
0.5027(58 
0.5630(53 
0.6483(49 
0.7219(62 
0.8319(39 
0.9439(21 
1.0039(21 
1.0044(30 



X4 

0.551( 38 
1.246(121 
1.662( 84 
2.199(127 
2.315(122 
1.762(154 
1.140( 98 
0.857( 29 
0.594( 19 
0.543( 12 
0.515( 7) 



Xi 



0.803(210 
0.915(225 
1.348(122 
1.786(196 
1.001( 94 
0.935( 80 
0.809( 87 
0.667( 75 
0.571( 56 
0.576( 36 
0.543( 22 
0.541( 12 



0.199( 5) 
0.561(22) 
0.722(25) 
0.968(39) 
1.189(48) 
1.117(65) 
1.104(58) 
1.013(25) 
0.754(16) 
0.631(13) 
0.536( 7) 



xi 



0.244(18) 
0.310(36) 
0.563(40) 
0.963(112) 
0.843(49) 
0.816(48) 
0.879(66) 
0.900(63) 
0.731(43) 
0.591(32) 
0.549(11) 
0.563(25) 



TABLE II: The data on quadratic and quartic fluctuations of 
light (it, d) and strange (s) quarks obtained from calculations 
on 16 3 x 4 (upper table) and 24 3 x 6 (lower table) lattices, 
respectively. 



dence than \2 an d xl separately. In fact, the difference 
between the N T = 4 and 6 results for X2/X2 i s we ^ un ~ 
derstood in terms of the small shift (cut-off dependence) 
of the transition temperature determined for these two 
different lattice sizes. 

We find that X2/X2 1S about 0.6 at T c and approaches 
unity at about 1.7T C . In both figures we show a band in- 
dicating the result for fluctuations in a non-interacting, 
massive Fermi gas. As can be seen this can explain the 
smaller fluctuations of strange quarks relative to light 
quarks only at temperatures larger than about 300 MeV, 
i.e. for T>1.5 T c . In quenched QCD and QCD with n f 
massless quark flavors the approach to the high tempera- 
ture limit has also been analyzed within HTL-resummed 
perturbation theory [2l| as well as in a straightforward 
perturbative calculation performed up to 0(g 6 \ng) [22T |. 
In both cases the results resemble closely properties of 
a Fermi gas with a temperature dependent mass term. 
Both approaches, however, use additional phenomeno- 
logical input to either re-sum next-to-leading order cor- 
rections to the HTL mass or fix the unknown scale for 
the 0(g 6 lng) correction. The lattice results presented in 
Fig. Q] overshoot the HTL result by about 5%. This is of 
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FIG. 1: The quadratic fluctuations of light and strange quark 
number versus temperature (upper figure) and the ratio of 
s and u quark fluctuations on lattices of size 16 3 x 4 and 
24 3 x 6, respectively. The solid lines show results for an ideal 
Fermi gas with mass m = 100 MeV and 200 MeV, respec- 
tively (top to bottom). In the lower figure we also show par- 
tially quenched results obtained with unimproved staggered 
fermions on 16 3 x 4 lattices These results have been given 
in Ref. [l2| as function of T/T c . This scale is shown on the 
upper x-axis in the lower part of the figure. We have used the 
mean value of the transition temperatures on the N T = 4 and 
6 lattices, i.e. To = 200 MeV, to compare with calculations 
performed with the p4-action. 



similar magnitude as the cut-off dependence seen in the 
current lattice results when comparing the N T = 4 and 
6 data in the temperature interval 1.5<T/T C <2. In view 
of the large change between the leading and the next to 
leading order HTL calculations as well as the systematics 
of the cut-off dependence seen in the lattice results, this 
seems to be a reasonably good agreement. 

We note that the results shown in Fig. [1] are in agree- 
ment with calculations performed with another improved 
staggered fermion action d, Q , the asqtad action. The 
results obtained with improved (p4 and asqtad) actions 
are, however, in contrast to calculations performed with 
the standard staggered action in quenched [llj and 2- 
fiavor [12| QCD. These calculations led to significantly 
larger values for \2 , which is understood in terms of the 
large cut-off effects in the standard discretization scheme 
(see appendix). In the absence of results on larger lat- 
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tices, which would have allowed for a continuum extrap- 
olation, the numerical results obtained within the stan- 
dard discretization scheme (IH [T3 | have been normalized 
to the corresponding ideal gas value evaluated also on lat- 
tices with finite temporal extent. It is well known that 
at temperatures a few times the transition temperature 
this procedure, which is correct in the limit of infinite 
temperature, over-estimates the actual cut-off distortion 
of numerical results. The normalized values thus end 
up to be substantially smaller than the results obtained 
with improved staggered fcrmion actions for which an ad- 
hoc normalization has not been performed. We discuss 
cut-off effects and their quark mass dependence for the 
standard and p4-action in somewhat more detail in the 
appendix. 

The situation is different for the ratio X2/X2 shown 
in the lower part of Fig. [TJ Here unknown normaliza- 
tion factors drop out. Nonetheless, also in this case re- 
sults obtained in calculations with improved actions are 
significantly larger than those obtained with a standard 
action. This difference may be due to the fact that in 
Ref. [T3| the strange quark sector has only been treated 
in the quenched approximation. 



IV. FLUCTUATION OF CONSERVED 
CHARGES 

In Fig. [5] we show results for quadratic (X2) and quar- 
tic {Xa ) fluctuations with X = B, Q and S. As can be 
seen, in all cases the quadratic fluctuations rise rapidly in 
the transition region where the quartic fluctuations show 
a maximum. This maximum is most pronounced for the 
baryon number fluctuations but is still visible also in fluc- 
tuations of the strangeness number. 

We note that results obtained on lattices with tem- 
poral extent N T = 6 are in good agreement with those 
obtained on the coarser N T = 4 lattice. The slight shift 
towards smaller temperatures, visible for the N T = 6 
data relative to the N T — 4 data, is consistent with find- 
ings for the equation of state, e.g. the trace anomaly 
(e — 3p)/T 4 , and also reflects the shift in the transition 
temperature observed when comparing the location of 
cusp in the chiral susceptibility [16]. The generic form 
of this temperature dependence, a smooth crossover for 
quadratic fluctuations and a peak in quartic fluctuations, 
is in fact expected to occur in the vicinity of the chiral 
phase transition of QCD. At vanishing chemical poten- 
tial and for vanishing quark mass the singular behavior 
of fluctuation observables like, e.g. the baryon number 
cumulants, is expected to be controlled by the singular 
part of the free energy that has the same universal struc- 
ture as that of a three dimensional 0(4) spin model 4 . 
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0.0301( 9) 
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0.0831(25) 
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0.111(25) 


0.599(41) 
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0.2694(23) 


0.6154(38) 
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0.2978( 5) 


0.6480(18) 


0.060( 5) 
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0.3313(18) 


0.6918(12) 


0.034( 6) 
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305.2 
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T[MeV] 


X? 


X? 


xf 


xi 


173.8 


0.0395(36) 


0.1517(30) 


0.094(56) 


0.215(10) 


179.6 


0.0581(48) 


0.1980(34) 


0.101(54) 


0.260(36) 


185.6 


0.0928(30) 


0.2700(47) 


0.102(32) 


0.435(34) 


195.0 


0.1600(37) 


0.3930(78) 


0.142(30) 


0.577(69) 


201.3 


0.2088(31) 


0.4928(46) 


0.047(18) 


0.346(29) 


204.6 


0.2299(41) 


0.5266(38) 


0.091(15) 


0.270(20) 


211.1 


0.2528(28) 


0.5723(32) 


0.067(22) 


0.237(18) 


224.3 


0.2794(21) 


0.6058(37) 


0.036(17) 


0.218(19) 


237.7 


0.3000(20) 


0.6388(15) 


0.043(23) 


0.156( 7) 


278.4 


0.3232(22) 


0.6637(16) 


0.029( 6) 


0.139( 7) 


362.9 


0.3354(11) 


0.6798(17) 


0.021( 7) 


0.135( 5) 


415.8 


0.3362(10) 


0.6777(12) 


0.023( 4) 


0.131( 2) 



TABLE III: The data on quadratic and quartic fluctuations 
obtained from calculations on 16 3 x 4 (upper table) and 24 3 x 6 
(lower table) lattices, respectively. 

This singular part of the free energy, f s , is controlled by 
a reduced 'temperature' t that is a function of tempera- 
ture as well as the quark chemical potentials \ix [H, [HJ . 
The latter add in even combinations to the reduced tem- 
perature in order to respect charge symmetry at /ix = 0. 
Baryon number cumulants are thus expected to scale like, 

m rri 2 — n — a 



with a ~ —0.25 [—0.015] [25[ denoting the critical ex- 
ponent characterizing the non-analytic structure of the 
specific heat in 3-d, 0(4) [0(2)] symmetric spin models. 
Starting at 6th order the baryon number cumulant thus 
will diverge at T c in the chiral limit. It has been argued 
in [7[ that this singularity will also show up in cumulants 
of the electric charge. 

On lattices with temporal extent N T = 4 and 6 and 
for the values of quark masses used here the transition 
temperature is close to 200 MeV. From Fig. [5] we thus 



4 As chiral symmetry is explicitly broken in numerical calculations 
with staggered fermions, the relevant symmetry group is, in fact, 



expected to be 0(2). 
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FIG. 2: Quadratic and quartic fluctuations of baryon num- 
ber, electric charge and strangeness. All quantities have been 
normalized to the corresponding free quark gas values, Xn' SB 
given in Table HVl 
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FIG. 3: The 6th order cumulant of electric charge and 
strangeness fluctuations evaluated on lattices of size 16 3 x 4. 

tions of the light u and d quarks. Only for temperatures 
T>1.7T C do the fluctuations of all three charges agree. 
This resembles the pattern of the relative strength of light 
and strange quark fluctuations seen in Fig. [T] 

On the coarser N T — 4 lattices we also have calculated 
the 6th order cumulants, x& ■> which vanish in the infi- 
nite temperature, ideal gas limit as well as in leading or- 
der perturbation theory [21( . For strangeness and charge 
fluctuations we show results in Fig. [3] Note that these 
6th order cumulants change sign at T ~ T c . Below T c 
they rise rapidly and reach a maximum at T ~ 0.95T C ; 
they approach zero from below and almost vanish for 
T>1.5T C . 



X 


B 


Q 


S 


X,SB 
X2 

xf SB 


1/3 
2/9tt 2 


2/3 
4/3tt 2 


1 

6/tt 2 



TABLE IV: Ideal gas values for quadratic (x2)and quartic 
(Xa) fluctuations of baryon number (B), electric charge (Q) 
and strangeness (S). 



conclude that at temperatures of about 1.5T C and larger 
quadratic and quartic cumulants of the fluctuations of 
B, Q and S are close to those of an ideal, massless quark 
gas, for which the pressure is given by 



SB 



P 



E 

f—u,d,s 



III 

60 



2 V T J 



1 

4^2 



.(8) 



Using the relations given in Eq. [T] one easily derives the 
corresponding ideal gas values for quadratic and quartic 
fluctuations of conserved charges. These are summarized 
in Table [TV] 

The upper part of Fig. [2] shows that over a wide 
temperature range quadratic fluctuations of strangeness 
are suppressed relative to those of baryon number and 
charge, which receive contributions also from fluctua- 



V. RATIOS OF CUMULANTS 

At low temperatures hadrons are the relevant degrees 
of freedom. The hadron resonance gas (HRG) model has 
been shown to provide a g ood description of thermal con- 
ditions at freeze-out [2g, |27|. Also fluctuations of the 
thermal medium have been successfully described in the 
framework of a HRG model 28]. Nonetheless, in a HRG 
model cumulants are monotonically rising functions of 
the temperature, while the 6th order cumulants change 
sign at T c . This indicates that a straightforward HRG 
model has to break down in the vicinity of the transition 
temperature. 

The partition function of the hadron resonance gas can 
be split into mesonic and baryonic contributions, 



V 



HRG 



l Ti = yjv E lnZ 7n t ( T ^ V ^B,^Q,^s) 



i(E mesons 



— ^2 In Z^.(T,V,fj, B ,IJ>Q,fJ>s) ,(9) 



y T 3 



£<G baryojis 



where 



lnZZ /B =T^ r*fc 2 ln( 
27r Jo 



1 T- z%e 



(10) 
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FIG. 4: The ratio of fourth and second order cumulants 
of baryon number (top), strangeness (middle) and electric 
charge (bottom) fluctuations. In the latter case we show 
curves for a HRG model calculated with physical pion masses 
(upper curve), pions of mass 220 MeV (middle) and infinitely 
heavy pions (lower curve). 

with energies e% = \Jk 2 + mf, degeneracy factors di and 
fugacities 

Zi = exp((Si/x B + Qip,Q + Sins)/T) . (11) 

The HRG model has been compared to lattice results 
previously [28] . In these earlier comparisons with lattice 
calculations, which had to be performed with rather large 
quark masses and without dynamical strange quarks the 
masses contributing to the HRG model were adjusted 
and strange contributions were suppressed. Here we no 



longer follow this strategy but compare directly the lat- 
tice results with an HRG model that is based on the ex- 
perimentally observed spectrum. We use the same HRG 
model as it also is used in heavy ion phenomenology and 
the analysis of freeze-out conditions [26|, l2jj • This ansatz 
for the HRG includes all meson and baryon masses from 
the particle data book with masses < 2.5 GeV that 
are characterized at least as '4 star states'. 

Also on the lattice with temporal extent N T = 6 cal- 
culations at low temperature are still performed on quite 
coarse lattices. We thus cannot expect to reproduce de- 
tails of the hadron spectrum well under these conditions. 
Nevertheless, we feel that it is now most appropriate 
to compare lattice results directly with the unmodified 
continuum version of the HRG model that is commonly 
used in the phenomenological treatment of QCD thermo- 
dynamics. By improving lattice calculations systemati- 
cally one will then be able to judge to what extent the 
HRG model does describe the thermodynamics of QCD 
in the hadronic phase. In fact, the current calculations 
of the QCD equation of state, e.g. the pressure and trace 
anomaly, (e — 3p)/T 4 , show deviations from the HRG at 
temperatures below the transition region. This is still 
the case in calculations on lattices with temporal extent 
N T = 8, i.e. closer to the continuum limit [2^, [3(|. We 
thus concentrate here on observables which are less sen- 
sitive to details of the hadron mass spectrum and em- 
phasize the charges of the relevant degrees of freedom 
contributing to the fluctuations. These are in general ra- 
tios of cumulants of Nb, Nq, and Ng. In fact, in the 
framework of a HRG model it is easy to convince oneself 
that the ratio of fourth and second order cumulants of 
baryon number fluctuations is completely independent of 
the actual value of hadron masses; 

all hadrons are heavy on the scale of the temperature. In 
that case the fugacity expansion of Eq. [T0l 

kz* = ^-(~) £(-i)*Hr a * a (wr) 

xcosh(^ B /T) , (12) 

is well approximated by its leading order term, i.e. the 
Boltzmann approximation. In this case the baryonic 
contribution to the pressure of a HRG is proportional 
to cosh(fi B /T) and the ratio xf/xf thus becomes in- 
dependent of the hadron mass spectrum. This is rea- 
sonably well reproduced by the lattice results shown in 
Fig. BJtop) ■ Note, however, that in the chiral limit it is 
expected that the cusp in \ f (Fig- EJ) is expected to be- 
come more pronounced and thus more prominent also in 
the ratio Xi 1X2- How strong this effect will be requires 
more detailed studies and also a better control over the 
continuum limit. Model calculations yield quite a differ- 
ent strength for the peak in xi / X2 • 

Even within the Boltzmann approximation the struc- 
ture of a HRG is, however, more complicated in the 
strange and/or electrically charged sectors. In these cases 
multiply strange hadrons or hadrons with charge Q = 2 
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contribute to the HRG. This enhances quartic fluctua- 
tions relative to quadratic ones and leads to a deviation 
of cumulant ratios from unity. This qualitative feature is 
indeed seen in the results obtained for xf I xf an d Xa Ixi 
shown in Fig. [U respectively. 

We note that the lightest hadrons, the pions, only con- 
tribute to cumulants of electric charge fluctuations. For 
these light states the Boltzmann approximation is not 
sufficient at temperatures close to T c . We thus have used 
the complete bosonic fugacity sum in the pion sector. 
These light hadronic states are also not well taken into 
account in current lattice calculations which still are be- 
ing performed with light quark masses that are about a 
factor two larger than the physical up and down quark 
masses. Moreover, the pion spectrum is distorted in cal- 
culations with staggered fermions due to cut-off effects 
that explicitly break flavor symmetry at non-vanishing 
lattice spacing (see footnote 3). In the case of electric 
charge fluctuations we therefore have analyzed in more 
detail the contribution of pions to the fluctuations in a 
HRG. In Fig. QJbottom) we show results for HRG model 
calculations with physical pion masses, with a pion mass 
of 220 MeV, which corresponds to the lightest pseudo- 
scalar state formed with the quark mass values currently 
used in our calculations, and without the pion sector, 
i.e. for infinitely heavy pions. Although the largest 
part of the electric charge fluctuations arises from heavier 
hadrons, the ratio of cumulants is sensitive to the pion 
sector. Already pion masses that are 50% larger than the 
physical value drastically reduce the contribution of the 
pion sector. This effect is even more significant in higher 
order cumulants. In the 6th order cumulant half of the 
fluctuations can be attributed to the light pion sector. 
This is obvious from Fig. [5] where we show the ratio of 
6th and 2nd order cumulants for electric charge fluctua- 
tions. In a HRG model the ratio of cumulants calculated 
without the light pion sector is a factor two smaller than 
that calculated with the physical pions included. In the 
low temperature regime the lattice calculations are con- 
sistent with the former. We thus conclude that lighter 
quark masses and calculations closer to the continuum 
limit are needed to correctly represent in lattice calcu- 
lations higher order cumulants that are sensitive to the 
light hadron sector. Irrespective of this we find, however, 
that the occurrence of maxima in x§ close but below T c 
signal the breakdown of the HRG model close to T c . The 
ratio Xe /X2 starts dropping below T c and is consistent 
with zero at T c . 

A similar behavior is found for higher order cumulants 
of strangeness fluctuations. We find that for both val- 
ues of the lattice cut-off the ratio xf I xf overshoots the 
HRG values in the transition region and this also holds 
true for the ratio xf /xf evaluated on the coarse 16 3 x 4 
lattices. Of course, this requires confirmation through 
calculations with lighter quark masses on lattices closer 
to the continuum limit. It may suggest that the contribu- 
tion of even heavier, experimentally not well-established 
multiple strange hadrons, which are not included in the 
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FIG. 5: The ratio of 6th and 2nd order cumulants of electric 
charge fluctuations evaluated on lattices of size 16 3 x 4. 

current version of the HRG model, is of importance in the 
transition region. In general we find, however, that the 
HRG model gives a fairly good description of cumulants 
of the fluctuations of conserved charges up to tempera- 
tures close to the transition temperature. 



VI. CORRELATIONS AMONG CONSERVED 
CHARGES 

The analysis of cumulant ratios presented in the pre- 
vious section suggests that for temperatures above 1.5T C 
fluctuations of baryon number, strangeness and electric 
charge agree well with the corresponding fluctuations in 
a non-interacting gas of light and strange quarks. In 
order to further test whether the relevant quasi-particle 
degrees of freedom indeed can be assigned to quarks, the 
analysis of correlations between different charges is quite 
instructive [hJ EH • Results for correlations between Nb , 
Ns and Nq are shown in Figs. [B] and [7] For completeness 
we also include the HRG-prediction. We note that the 
model is not expected and, in fact, does not capture the 
fluctuations and correlations in the transition region cor- 
rectly. However, as shown in Figs. 6 and 7 as well as in 
Fig. 4 it reproduces qualitatively the features of fluctu- 
ations and correlations. In particular, the drop and rise 
seen in our numerical results for BQ and BS correlations, 
respectively, is seen also in the HRG model calculations. 
To put these observations on a more quantitative basis 
will require further calculations at lower temperatures. 

In Fig. [S] the correlation functions in the numerator 
project only on the charged baryon sector of the spec- 
trum. The ratio Xii lx-2 therefore approaches 5 1/2 in 
the low temperature limit as the numerator receives con- 
tributions from protons and anti-protons only, while the 



We ignore here small mass differences of charged and neutral 
hadrons. 
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T[MeV] 


BQ 
Xll 


Xii 


QS 

All 


176.0 


0.00630(35) 


-0.0175( 9) 


0.0374( 2) 


186.4 


0.01332(40) 


-0.0564(14) 


0.0811( 5) 


190.8 


0.01670(26) 


-0.0727(13) 


0.0966( 9) 


195.4 


0.02023(43) 


-0.0924(19) 


0.1133(14) 


202.4 


0.02811(32) 


-0.1369(14) 


0.1548(13) 


204.8 


0.02959(41) 


-0.1544(20) 


0.1687(15) 


209.6 


0.03024(58) 


-0.2089(19) 


0.2244(19) 


218.2 


0.02912(17) 


-0.2395( 8) 


0.2488( 9) 


253.7 


0.01625(12) 


-0.2988(16) 


0.3045( 6) 


305.2 


0.00765( 6) 


-0.3245(11) 


0.3303( 5) 


444.4 


0.00272( 1) 


-0.3289( 4) 


0.3312( 3) 


T[MeV] 


BQ 
All 


V s5 
All 


QS 
All 


173.8 


0.00830(88) 


-0.0228(23) 


0.0469( 8) 


179.6 


0.01040(120) 


-0.0373(26) 


0.0635(14) 


185.6 


0.01611(74) 


-0.0606(18) 


0.0868(18) 


195.0 


0.02406(49) 


-0.1119(33) 


0.1312(31) 


201.3 


0.02697(72) 


-0.1549(24) 


0.1739(18) 


204.6 


0.02641(72) 


-0.1771(31) 


0.1930(13) 


211.1 


0.02312(67) 


-0.2066(25) 


0.2208(16) 


224.3 


0.02309(90) 


-0.2332(25) 


0.2444(20) 


237.7 


0.01529(66) 


-0.2695(23) 


0.2812(10) 


278.4 


0.00529(77) 


-0.3126(10) 


0.3157( 8) 


362.9 


0.00145(63) 


-0.3325( 8) 


0.3357(10) 


415.8 


0.00122(34) 


-0.3338(13) 


0.3353(12) 



TABLE V: The data on correlations between different quan- 
tum number fluctuations obtained from calculations on lat- 
tices with temporal extent N T = 4 (upper table) and 6 (lower 
table), respectively. 
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FIG. 6: Correlations of electric charge and strangeness 
with baryon number normalized to quadratic fluctuations of 
baryon number. 



VII. CONCLUSIONS 



denominator also receives contributions from the neu- 
trons. The ratio X11 /x% 1 however, will approach zero in 
the low temperature limit as the lightest baryons carry 
no strangeness! 

In Fig. [7] we show the correlation among Q and S nor- 
malized to the quadratic fluctuations of electric charge. 
A similar ratio, where strangeness fluctuations have been 
used to normalize the Q-S correlations, has been dis- 
cussed in [l(| and has also been calculated in 2-flavor 
QCD with a quenched strange quark sector fl3| . We pre- 
fer the above normalization, as it emphasizes the non 
trivial correlations between Q and S that persist in the 
high temperature phase of QCD. Unlike xf the charge 
fluctuations X2 approach the ideal gas value rapidly 
above T c (see Fig. [5]). The deviations from ideal gas 
behavior seen in Fig. [7] thus are mainly due to the de- 
viations of x?i from ideal gas behavior. This is shown 
in the lower part of Fig. [Jj Apparently, the temperature 
dependence of x?i is very similar to that of x§ ■ 



We have analyzed the fluctuations of baryon num- 
ber, electric charge and strangeness in finite temperature 
QCD at vanishing chemical potential. A comparison of 
calculations performed with 0(a 2 ) improved staggered 
fcrmions for two different values of the lattice cut-off 
shows that these effects are generally small; cut-off ef- 
fects as well as the dependence of observables on the light 
quark masses become, however, more severe for higher 
order cumulants. 

We find fluctuations and correlations of conserved 
charges are well described by an ideal, massless quark 
gas already for temperatures of about (1.5 — 1.7) times 
the transition temperature. Deviations from ideal gas be- 
havior seem to be mainly induced by the strange quark 
sector for which the quadratic quark number fluctuations 
approach the ideal gas limit more slowly. This cannot 
only be explained by the quark mass dependence of an 
ideal gas but seems to reflect a significant quark mass 
dependence of the interaction at high temperature. 

At low temperature we find that fluctuations and cor- 
relations of conserved charges are reasonably well de- 
scribed by a hadron resonance gas up to temperatures 
close to the transition temperature. Higher order cumu- 
lants, however, signal that the resonance gas prescription 
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APPENDIX: QUARK MASS DEPENDENCE OF 
CUT-OFF EFFECTS 



Cut-off effects arising from the introduction of a finite 
lattice spacing and, in particular, their dependence on the 
discretization scheme have been extensively discussed in 
the ideal gas limit. However, standard as well as 0(a 2 ) 
improved staggered fermion formulations have usually 
been analyzed in the massless limit, which is appropriate 
for the light quark sector of QCD. For the strange quark 
sector non-zero mass effects may play a role, in partic- 
ular close to the QCD transition temperature which is 
only about twice as large as the strange quark mass. We 
present here an analysis of the quark mass dependence 
of cut-off effects in the ideal gas limit. 

The quark number susceptibilities are obtained as the 
second derivative of the pressure, p/T 4 , with respect to 
the chemical potential normalized to the temperature, 
fi/T. An analysis of the cut-off and quark mass depen- 
dence of X2 ,S t nus can follow closely the corresponding 
analysis performed for the pressure 6 in the massless limit 
32]. On a lattice of size N^N T this is given by 



has to break down at temperature values close but below 
T c . 

The current analysis has been performed with light 
quarks that are one tenth of the strange quark mass. We 
have shown, that higher order cumulant ratios like, for 
instance x& /X2' become quite sensitive to the pseudo- 
scalar Goldstone mass. In numerical calculations with 
staggered fermions this also means that results become 
sensitive to a correct representation of the entire Gold- 
stone multiplet. Calculations with smaller quark masses 
closer to the continuum limit will thus be needed in the 
future to correctly resolve these higher order cumulants, 
which will give deeper insight into the range of appli- 
cability of the resonance gas model at low temperature 
and the non-perturbative features of the QGP above but 
close to T c . 
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|i = 2 (^) 3 £ ln( ^ + m2) ' ( A -D 

^ a ' p,p4 

with, D 2 = £* =1 Dl where 

Di = c w sin(pi) + c 30 sin(3pi) + 2c i2 sin(pi) x 

[cos(2p 2 ) + cos(2p 3 ) + cos(2p 4 )] (A.2) 

and D% , D3 and D4 are obtained by cyclic permutation of 
the piS. For the standard (naive) discretization scheme, 
one has do = U C30 — C12 = 0; the discretization scheme 
used in the p4-action corresponds to do = 3/4, c 30 = 0, 
d2 = 1/24 and for the Naik action it is do = 9/8 and 
C30 = —1/24 and d2 = 0. 

We evaluate here quark number susceptibilities as a 
function of quark mass on lattices with finite tempo- 
ral extent N T and infinite spatial volume for the stan- 
dard (naive) staggered fermion discretization and the im- 
proved (p4) discretization scheme. For vanishing quark 
mass results on the cut-off dependence of quark number 



6 Note that at vanishing chemical potential the pressure is simply 
related to the free energy density, p = — f. 
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FIG. 8: Cut-off dependence of quark number fluctuations in 
the infinite temperature, ideal gas limit. Shown are results 
for staggered fermions in the standard (naive), p4 and Naik 
discretization scheme. 



susceptibilities in the ideal gas limit have also been pre- 
sented in Ref. [H, [34]]. There also some results for X4 
and X6 can be found. 

Starting with the expressions for p/T 4 at non-zero 
chemical potential (p4 — > p4 — i/j,) and taking two deriva- 
tives with respect to fi/T one finds in the standard 
scheme for /i = 0, 



Xa(m, N T ) 

A 
B 



^(A-B + AB), 



2, 



2 sin (pn) 



m2 + Z)i=x sin 2 (p 4 )' 
2 cos 2 (p 4 ) 

m2 + Si=i sm2 (ft) 



(A.3) 



the corresponding expressions for the p4 and Naik ac- 
tions are somewhat more complicated but can be derived 
straightforwardly. 

In the continuum limit the quark number susceptibility 
is given by, 



T 3 7T 2 



dfc- 



X 2 (to/T,oo) 

with E 2 = k' 2 + 
quark mass X2 is quadratic in the mass, 



-E/T 



,2 ' 



(A.4) 



to (l + c- E / T ) 
i 2 . We note that for small values of the 



X2 (m/T, oo) = 1 - JL (^) 2 + 0((m/T) 4 ) . (A.5) 



We have evaluated X2( m , -^r) f° r various values of _/V r 
and fixed quark mass as well as for fixed N T varying the 
quark mass. The results for the naive and p4 actions are 
shown in Fig. [5] 

As expected, the cut-off dependence of %2 at non- 
zero values of the quark mass closely follow the pattern 
known from the analysis of the pressure at finite N T and 
m/T = 0. While the standard discretization scheme 
shows large deviations from the continuum result even 
for N T ~ 10, the results for the p4 action come close to 
the continuum result already for N T = 6. This also is 
true for the ratio X2{m, N T ) /x2(0, N T ). Although in this 
ratio a large part of the cut-off dependence cancels, the 
quark mass dependent corrections remain. 

We note, that for finite N T the quark mass depen- 
dence of X2("ij At)/x2(0, N t ) in the naive discretization 
scheme is weaker than in the continuum limit. This is 
qualitatively different from the differences seen in the 
lower part of Fig. [T] between the standard and improved 
discretization schemes. 
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